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ABSTRACT 

We present the first direct measurement of the mean Halo Occupation Distribution (HOD) of X-ray 
selected AGN in the COSMOS field at z < 1, based on the association of 41 XMM and 17 C- 
COSMOS AGN with member galaxies of 189 X-ray detected galaxy groups from XMM and Chandra 
data. We model the mean AGN occupation in the halo mass range logM2oo[MQ] = 13 — 14.5 with 
a rolling-off power-law with the best fit index a — 0.06(— 0.22; 0.36) and normalization parameter 
f a = 0.05(0.04; 0.06). We find the mean HOD of AGN among central galaxies to be modelled by a 
softened step function at \ogM h >logM mi „ = 12.75(12.10, 12.95)M Q while for the satellite AGN HOD 
we find a preference for an increasing AGN fraction with M/j suggesting that the average number of 
AGN in satellite galaxies grows slower (a s < 0.6) than the linear proportion (a s = 1) observed for the 
satellite HOD of samples of galaxies. We present an estimate of the projected auto correlation function 
(ACF) of galaxy groups over the range of r p — 0.1 - 40 h _1 Mpc at (z) = 0.5. We use the large-scale 
clustering signal to verify the agreement between the group bias estimated by using the observed 
galaxy groups ACF and the value derived from the group mass estimates. We perform a measurement 
of the projected AGN-galaxy group cross-correlation function, excluding from the analysis AGN that 
are within galaxy groups and we model the 2-halo term of the clustering signal with the mean AGN 
HOD based on our results. 

Subject headings: Surveys - Galaxies: active, clusters - X-rays: general - Cosmology: Large-scale 
structure of Universe - Dark Matter 



1. INTRODUCTION 

In the past decade it has become clear that the major- 
ity of galaxies have likely had one or more brief active 
periods, shifting our view of astrophysical black holes 
(BHs) from the role of exotic phenomena to fundamental 
ingredients of cosmic structures. The question of when 
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and under which physical conditions the active galactic 
nuclei (AGN) activate is important for understanding not 
only the origin and evolution of BHs but also the origin 
and evolution of galaxies. 

The clustering analysis can powerfully test theoretical 
model predictions and address which physical processes 
trigger AGN activity. In the framework of the cold 
dark matter structure formation (CDM), the spatial 
distribution of AGN in the Universe (described by the 
correlation function, CF) can be related to how AGN are 
biased with respect to the underlying matter distribution 
and to the typical mass of dark matter halos (DMHs) 
in which they reside. There have been several studies of 
the bias evolution of optical quasars with redshift, based 
on survey such as 2QZ and SPSS (ICroom et all [ 2005; 
Porciani 'fc Norberd [200l IShen et all 120091: IRoss et al.l 
2009). All the previous studies infer the picture that the 
quasar bias evolves with redshift following a constant 
mass evolution in the range logM^h^M©] ~ 12.5 — 13, 
i.e. halo masses similar to group scales, where the com- 
bination of low velocity dispersion and moderate galaxy 
space densi ty yields to the highest probability of a clos e 
encounter (jHopkins et al.l 120081: IMcIntosh etaLI [2009) . 
Models of major mergers between gas-rich galaxies 
appear to naturally produce many observed properties 
of quasars, such as the shape and the evolution of the 
quasar luminosity function and the large-scale quasar 
clustering as a fun c tion o f luminosity and redshift 
(Hop kins et all [2007L [200l IShenl l2lM IShankar et all 
120091 120101: iBonoli et al.l 12009ft . supporting the sce- 
nario in which major mergers dominate the bright 
quasar populations. On the other hand, the major- 
ity of the results on the clustering of X-ray selected 
AGN, suggest a picture where moderate-luminosity 
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AGN live in massive D MHs (12.5 < logM h \h.- 1 M, 



< 13 5) up to z ~ 2 (iGilli et all 120051: lYang et al 
2001 IGilli et all l200l IHickox et al.1 1200a ICoil et al 



2003: iKrumpe et all I2010L ICappelluti et al.1 120101: 



Krumpe et alj|2010l; (Allevato et al.1 120111: IMivaii et alj 
201 ID . i.e. X-ray selected AGN samples appear to 



cluster more strongly than luminous quasars (sec 
ICappelluti. Allevato fc Finoguenovl [2012, for a review 
on the argument). The reason for this is not completely 
clear, but several studies argued that this large bias and 
DMH masses could suggest a different AGN triggering 
mechanism with respect to bright quasars characterized 
by galaxy merger-induced fueling. Several studies on 
the morphology of the AGN host galaxies have demon- 
strated that major mergers of galaxies are not likely to 
be the single dominant mechanism r esponsible for trig- 
gerin g AGN activity at low ( z ~ I) dGeorgakakis et al.1 



20071: ISilverman et al l [20091: lG eorgakakis et al.1 l200£ 
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Dunlop et all 120031: iGrogin et all 2005; Pierce et al 
20071: [Gabor et all 120091: iReichard et al.ll2009f iTal et al 
2009t iCisternas et al.| 120111: ISilverman et al.M201lD and 
high redshift (z ~ 2) (iRosario et al.ll2011l : lKocevski et al.l 
120111: iSchawinski et al.ll2011[ ). 

The theoretical understanding of galaxy clustering and 
bias factor has been greatly enhanced through the Halo 
Occupation Distribution (HOD) framework. In this 
framework, the virialised DMH with typical overdensi- 
ties of A 2 oo (defined w.r.t. the mean density) are de- 
scribed in terms of the probability P(N\AI) of a halo 
of given mass M of having N galaxies. A simple way to 
model the complicated shape of N(M) is by assuming the 
existence of two separate galaxy populations within ha- 
los, central and satellite galaxies. This m ethod has been 
used extensively to interpret galaxy CFs (lHamana et al. | 
2004 [Tinker et al.ll2005l lPhleps et aTll2006l: iZh eng et al.l 
2007t IZehavi et al.ll20ia Izheng et al.ll2009D to constrain 
how various galaxy samples are distributed among DMH 
as well as whether these galaxies occupy the centers of 
the DMHs or are satellite galaxies (iKravtsov et al.ll2004J: 
Zhen get al.1 120051 : IZehavi et al.1 120051 : iRichardson et al.1 
20121 ). These two populations can be modeled with an 
HOD described by a step function above a halo mass 
limit for central galaxies, and a power-law for satellite 
galaxies. 

Similarly, the problem of discussing the abundance and 
spatial distribution of AGN can be reduced to studying 
how they populate their host halos. In fact the observed 
departure of the AGN CF from a power law on small 
scales (1-2 h^Mpc) can be physically interpreted in the 
language of the halo model, as the transition between two 
scales - from small scales lying within the DMHs (1-halo 
term) to those larger than the halo (2-halo term). The 
1-halo term constrains the HOD of satellite AGN and 
gives us the average profile of pairs of AGN in groups 
and clusters of galaxies. The 2-halo term reflects the 
large-scale AGN bias driven by the typical mass of the 
hosting halos. 

Due to the low number density of AGN, there have 
been few results in the literature studying the AGN cor- 
re lation function using: HOD modeling. Pre yious works 
of iPadmanabhan et al.l ([20091) at z < 0.6 and lShen et al.1 
([2010D at z = 3 - 4 on QSO using the HOD modeling 
found that > 25% and ^10% of th eir QSOs, respectively, 
are satellites. IMivaii et al.l ([201 1[ ) described for the first 



time the shape of the HOD of X-ray selected AGN. By us- 
ing the cross correlation function of ROSAT-RASS AGN 
with SDSS galaxies, they modelled the mean AGN oc- 
cupation of DMHs suggesting that the satellite AGN 
fraction increases slow (or may even decrease) with Mh, 
in contrast with the satellite HOD of luminosity-limited 
samples of galaxies. Cosmologi cal hydrodynamic simu- 
lations have been performed in iChatteriee et al.l ([20121 ) 
to study the mean occupation function of low-luminosity 
AGN as function of redshift and luminosity. They used 
a softened step function for the central component plus 
a rolling-off power-law for the satellite component with 
a = 0.3 — 1.4 depending on the redshift and AGN lu- 
minosity. Their results suggest a strong evolution of the 
AGN occupancy in the redshift range z — 1 — 3 estimated 
at three different luminosities L B ol > 10 38 , 10 40 , 10 42 
s _1 erg. IRichardson et al.l ([2012D modelled the HOD of 
SDSS quasars at z ~ 1.4 following this parametrization. 
They found that the satellite occupation becomes signif- 
icant at mass ~ 10 14 h _1 MQ, i.e. only the most massive 
halos host multiple quasars at this redshift and only a 
small satellite fraction (f sat = 7.4 ± 1.3 x 10~ 4 ) of SDSS 
quasars is required to fit the clustering signal at small 
scales. Moreover, they measured that the quasars HOD 
steepens considerably going from z=1.4 to 3.2 over halo 
mass scales 10 13-14 h _1 M Q and that the characteristic 
halo mass increases with z for central quasars. 

Despite the diverse methods for studying the HOD, 
counting the number of AGN within galaxy groups can 
constrain quite directly the average AGN number within 
a halo as a function of halo mass. The total mass of 
galaxy groups can be estimated via gravitational lensing 
and the distribution of AGN within halos can be investi- 
gated in groups by means of the distribution of the AGN 
host galaxies. Separating the contribution to the occu- 
pation of halos from AGN in satellite or central galax- 
ies can advance our understanding of the co-evolution 
AGN/galaxy and is related to the mechanism of AGN 
activation. 

On the other hand, the cross-correlation (CCF) of 
AGN with galaxy groups provides additional information 
about how galaxies and BH co-evolve in dense environ- 
ments. In fact the physical processes that drive galaxy 
evolution, such as the available cold gas to fuel star for- 
mation and the BH growth, are substantially different in 
groups and clusters compared to the field. Many stud- 
ies over the past decade have presented an evidence that 
AGN at z ~ 1 are more frequently found in groups com- 
pared to galaxies (Georgakakis et al. 2008, Arnold et 
al. 2009). X-ray observations reveal that a significant 
fraction of high-z clusters show overdensities of AGNs in 
their outskirts (Henry et al. 1991, Cappi et al. 2001, 
Ruderman et al. 2005, Cappelluti et al. 2005). 

In this paper, we perform the first direct measure- 
ment of the mean halo occupation of X-ray AGN and 
of the projected cross-correlation function of AGN with 
galaxy groups, using a sample of X-ray selected AGN 
and galaxy groups in the COSMOS field at z < 1, from 
XMM and Chandra data. We use a ACDM cosmology 
with Q M = 0.28, fl A = 0.72, Q b = 0.045, a s = 0.8. 
For comparison with previous measurements we refer to 
correlation lengths and distances in units of Mpc 
comoving, where Ho = 100 h km s _1 Mpc -1 . AGN lu- 
minosities and galaxy groups masses are calculated using 
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Fig. 1. — Left panel: X-ray rest-frame soft luminosity Lx as a function of redshift for XMM-AGN (circles) and C-COSMOS selected 
AGN (triangle) with know spec-z (filled) and or phot-z (open). The open green squares indicate XMM-AGN within galaxy groups, while 
the open magenta hexagons represent 17 additional C-COSMOS sleeted AGN in galaxy groups. Right panel: Galaxy group masses Af200 
as a function of redshift for the whole galaxy group sample (black circles) while open circles (squares) indicate XMM-AGN (C-COSMOS 
AGN) hosting groups. The mass estimates are defined with respect to 200 times the mean density, in units of Mq. 



h = 0.72. 

2. THE CATALOGS 
2.1. AGN Catalog 

The Cosmic Evolution Survey (COSMOS) is a mul- 
tiwavelength survey over 1.4 x 1.4 deg 2 of equa- 
torial sky observed by the most advance d astro- 
nomica l facilities HST (IScoville et all [20071) . SUB- 
ARU dTaniguchi et all 120071), Spitzer ^Sanders et all 
120071). GALEX (IZamoiski et all 120071) XMM - Newto n 
(IHasinger etHI 120071; ICappelluti eTal! [20071 120091) . 
Chandra l|Elvis et all 12009; IPuccetti et all 120091) . de- 
signed to thoroughly probe the evolution of galaxies, 
AGN, and dark matter in the context of their cosmic 
envir onment (LSS). XMM-Newton surveyed 2.13 deg 2 of 
this field in the 0.5-10 keV energy band for a total of ~ 
1.55 Ms providing a large sampl e of 1822 point-like X-ray 
sources (jCappelluti et al.ll2009l ). down to limiting fluxes 
of ~ 5 x 10" 16 , ~ 3 x 10" 15 and ~ 7 x 10" 15 erg cm^s- 1 
in the soft, hard and ultra hard band, respectively. 

The inner part of the COSMOS field (~ 0.92 deg 2 ) 
has been imaged for a total of 1.8 Ms by Chandra 
down to limiting fluxes 3-4 times deeper than XMM- 
COSMOS. Of the 1822 XMM sources, 945 have been ob- 
served by Chandra and 875 of them ar e presented in the 
C-COSMOS point-li k e source catalog (lElvis et al.ll2009t 
IPuccetti et all l2009t iCivano et all I2012I ). iBrusa et all 
42O1O0 presented the XMM-COSMOS multiwavelength 
catalog of X-ray sources with optical/near- infrared iden- 
tification, multiwavelength properties and redshift infor- 
mation. Starting from this catalog, we restricted the 
analysis to a sample of X-ray AGN (we removed nor- 
mal galaxies and ambiguous sources) detected in the soft 
band which guarantees the largest sample of X-ray AGN 
in the COSMOS field, compared to AGN observed in the 
hard or ultra-hard band. Moreover, the selection in only 
one band allows a more simple treatment of the AGN X- 



ray luminosity function used to correct the AGN HOD 
(© and of the XMM-and C-COSMOS sensitivity maps 
used to generate the AGN random catalog f £|4.3[) . 

Specifically, we selected a sample of 280 and 83 soft 
XMM-COSM OS AGN with spect roscopic and photomet- 
ric redshift (jSalvato et all 120111 ) z < 1, respectively. 
The use of photometric redshifts for group member- 
ship assignm e nt ha s been successfully demonstrated in 
iGeorge et all (|201lD . Note that 184/363 sources are also 
Chandra detected AGN. In order to test if the AGN halo 
occupation significantly changes including sources from 
the C-COSMOS catalog that are only Chandra detected 
(hereafter C-COSMOS AGN) (see jgjl, we included in 
the analysis a sample of 107 and 61 AGN, detected in 
the soft band, with known spectroscopic or photometric 
redshifts z < 1, respectively. The rest-frame soft X-ray 
luminosity as a function of redshift is shown in Fig. Q] 
(Left Panel) for XMM-COSMOS AGN (circles) and C- 
COSMOS AGN (triangles). 

2.2. Galaxy Group Catalog 

We used a catalog with 270 galaxy groups, de- 
tected in the co-added XMM-Newton and Chandra data. 
The general dat a redu ction process is described in 
iFinoguenov et all (|2007l) and details regarding improve- 
ments and modifi c ation to th e initial catalog are g iven in 
ILeauthaud eTail (pOlOl) and IGeorge et all (pOllh . The 
identification of the groups has been done using the 
red sequence technique, and spectroscopic identifica- 
tion of groups has been achieved t hrough zCOSMOS- 
B RIGHT program ([Lilly et alj|20090 . targ eted follow-up 
using IMACS/Magellan a nd FORS2/VLT (IGeorge et all 
l20Tll) . Gemini GMOS-S (|Balogh et all 1201 11) as well as 
through secondary targets on Keck runs by COSMOS 
collaboration. At z < 1 the spectroscopic completeness 
has achieved 90% (A. Finoguenov 2 12, in preparation). 

As described in ILeauthaud et all (|2010| ). the total X- 
ray fluxes have been obtained from the measured fluxes 
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by assuming a beta profile and by removing the flux 
that is due to embedded AGN point sources. Some 
of the faint Chandra AGN could not be removed from 
the group fluxes, with their contribution to the to- 
tal flux being < 10%. The rest-frame l umino sities 
have been comp u ted i n iFinoguenov et al.l (|2007| ) and 
ILeauthaud et al.l ((2010D from the total flux following 
£o.i-2.4fc e v = 4Trd 2 L K{z,T)Cp(z,T)Fd, where K(z,T) is 
the K-correction and Cp(z,T) is an iterative correction 
factor, while to estimate the tem perature of e ach g roup 
they used the Lx — T relation of iMarkevitchl (|1998f ). 

A quality flag (hereafter 'XFLAG') is assigned to the 
reliability of the optical counterpart, with flags 1 and 
2 indicating a secure association, and higher flags indi- 
cating potential problems due to projections with other 
sources or bad photometry due to bright stars in the fore- 
ground. In detail, XFLAG=1,2 are assigned to groups 
with a confident spectroscopic association, while systems 
wit h only the red sequence i dentification have XFLA G=3 
(see ILeauthaud et al.H2010L for more details). 

The line-of-sight position of the group is assigned to 
be the centroid of the X-ray emission (the accuracy 
of the determination of the X-ray center is higher for 
XFLA G—l ) . If the X-ray centroid is not precise enough 
to be used directlythe Most Massive Central Galaxy 
(MMCG) located near the peak of the X-ray emission 
has been us ed to trace the center of the DM halos of 
groups (see ILeauthaud et al.l 120101: iGeorge et all 120111 
for more details). Group masses M200 are assigned 
from an empirica l mass -luminosity relation, described in 
ILeauthaud et al.l (|2010D . 

log w (M 2 <n,c) =Pq- log w E(z) + log w (M ) (1) 
+Pi[log w (L x /E(z)) - log w (Lo)} 

where M200 is the mass within the radius containing the 
density of matter 200 times the critical density, in units 
of M . {p ,pi} = {0.729538,0.561657} are the fitting 
parameters, {logMo, logLo} = {13,42.5} are the calibra- 
tion parameters and E(z) is the correction for redshift 
evolution of scaling rel ations, which has been shown in 
ILeauthaud et al.l ((2010) to reproduce well the Lx — M 
relation of COSMOS groups. 

In order to be consistent in comparing these mass val- 
ues with the ones obtained studying the clustering prop- 
erties of groups, we accounted for the difference between 
the mass defined with respect to 200 times the crit- 
ical density and with respect to 200 times the mean 
density (hereafter M 2 oo refers to masses obtained us- 
ing the definition with respect to the mean density). 
In fact, the absolute bias which we are going to derive 
from the DM correlation function is based on the shape 
of the DM mass function defined with respect to the 
mean density. Starting from the Navarro-Frenk- White 
(NFW) profile with a concentration parameter c = 5, we 
derived the relation between the two mass definitions, 
M 200 ,m = M 200 ,c x n(z)-°- 13i . 

In this work, we make use of galaxy groups with z < 1 
and identification flag XFLAG < 3 which removes prob- 
lematic identification cases, obtaining a catalog of 189 X- 
ray galaxy groups over 1.64 deg 2 with a rest-frame 0.1-2.4 
keV luminosity range of 41.3 < log(Lx[s~ 1 erg}) < 44.1, 
and mass range of 13 < logM2oo[Mo\ < 14.5 (see Fig. [TJ 
Right Panel). 



2.3. AGN in galaxy groups 

We define group members as AGN located within < 
3(7 and < -R200 from the group centers, where a is the 
group line-of-sight velocity dispersion and -R200 is the 
virial radius of a group within which the mean density is 
200 times the mean density of the Universe at the group 
redshift. 

In this analysis we used the sample of 363 XMM- 
COSMOS AGN with z < 1 described in SJ2ZD and we 
found 41 sources (35/41 are also Chandra detected) 
in galaxy groups with median (z) — 0.55 and median 
(Lx) = lO 42 6 s _1 erg. When we include in the analysis 
the sample of 168 C-COSMOS AGN f ^TT]) . we found 17 
additional AGN in galaxy groups, with known spectro- 
scopic or photometric z < 1. As expected C-COSMOS 
AGN have lower soft fluxes respect to XMM-COSMOS 
AGN at any redshift, with a median Lx = 10 417 s _1 erg. 
In particular we found 2 galaxy groups at z ~ 0.1 and 
2 groups at z ~ 0.7 — 0.8 with 2 AGN per halo and 

1 group with 3 AGN at z ~ 0.35. All the properties 
of AGN and galaxy group samples are summarized in 
Table [Q while Table M shows the catalog of 58 XMM 
+ C-COSMOS AGN in galaxy groups. For the sources 
with known spectroscopic redshifts (47/58) we know the 
classification in B L or non-BL AGN as d escribed in 
iBrusa et al.l (|2010f ) and iCivano et ail (|2012f ) for XMM 
and C-COSMOS AGN, respectively. In detail, we found 
43 non-BL and 4 BL in the XMM + C-COSMOS AGN 
sample. Fig. Q] (Left Panel) shows the rest-frame soft X- 
ray luminosity as a function of redshift for the subsam- 
ples of XMM (green squares) and C-COSMOS selected 
AGN (magenta hexagon) within i?200- 

We cross-matched the sample of AGN in groups with 
a galaxy members hip catalog (|Leauthaud et al.l 120071 : 
IGeorge et al.|[20TT[) and we verified that AGN classified 
as group members based on our method, have host galax- 
ies associated with the same galaxy groups. We divided 
the sample of 58 AGN in groups in two subsets, accord- 
ing to their association with BCGs. In detail we found 
that 22/58 (16/41) AGN are in central galaxies, while 
36/58 (25/41) are in satellites. 

3. HALO MASS FUNCTION AND AGN HOD 

Fig. [2] (Left Panel) shows the mass function of all X- 
ray galaxy groups and those marked by AGN presence, 
showing separately the contributions of groups hosting 
an AGN in central or satellite galaxies. We calculated 
the mass function by using the standard 1/V max method 
(Schmidt 1968) and we counted twice galaxy groups with 

2 AGN. Hence, in the i th mass bin, the comoving space 
density (n{) and its corresponding error (<jj) are com- 
puted by (see Bondi et al. 2008): 




In estimating the average number of AGN occupying a 
halo of mass M200, some major effects need to be taken 
into consideration. The sample of AGN in -R200 is a flux- 
limited sample and brighter AGN are detected at higher 
redshift. Similarly, we observe galaxy groups with small 
halo mass only at low redshift. However the relatively 
small number of COSMOS AGN in galaxy groups does 
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Fig. 2. — Left Panel: Mass function of X-ray galaxy groups (triangles) and AGN host groups (circles). The black squares (red crosses) 
show the mass function of groups hosting an AGN in satellite (central) galaxies. Right Panel: Same as left panel when correcting for the 
AGN soft XLF and the redshift evolution. 



TABLE 1 

Properties of the Group and AGN Samples 



(1) 


(2) 


(3) 


(4) 


Sample 


N 


(,)• 


(Lx) h 


XMM-AGN 


363 


0.66 


10 42.8 


C-COSMOS AGN 


168 


0.56 


10 41.8 


XMM-AGN in R 20 o 


41 


0.55 


10 42.6 


C-COSMOS AGN in R 200 


17 


0.53 


1Q 41.7 


XMM-AGN in the field d 


253 


0.67 c 


10 42.8 


XMM+C-COSMOS AGN in R 200 


58 


0.55 


10 42.3 


- Satellites 


36 


0.56 


10 42.2 


- Centrals 


22 


0.50 


10 42.4 








(A/200 > c 


Galaxy Groups 


189 


0.56 


1Q13.60 


XMM+C-COSMOS AGN host groups 


52 


0.55 


^q13.62 



a Median redshift of the sample. 

In units of hj erg s — 1 
c Mass denned respect to 200 times the mean density known with a 
20% error, in units of Mq . 
^ AGN sample used in estimating the CCF. 
G Only spectroscopic rcdshifts. 

not allow us to select a volume complete subsample by 
using a cut in luminosity at logLx[erg s _1 ]=42.4. We 
have therefore made a correction for the effect of changes 
in the AGN density as a function of redshift and limiting 
luminosity, by using the AGN X-ray luminosity func- 
tion (XLF). It has been shown in different works that 
the Luminosity Dependent Density Evolution (LDDE) 
model for the XLF provides the best framework that 
describes the evolutionary properties of AGN, bo th in 
the soft (IMivaii et all 120001: lHasinger et al.l 120051) an d 
hard X-ravs (jUeda et al.l l2003t lLa Franca et al.l 120051) . 
In this paper we modelle d the AGN so ft XLF with the 
LDDE XLF described in lEbrero et~aT1 (|2009D . because 
they modelled the soft XLF including in the analysis both 



type 1 and type 2 AGN, but we verified that using differ- 
ent best-fit parameters of the global XLF expression, the 
resulting mean AGN occupation stays within the error 
bars. 

Then for each AGN redshift we defined two weights: w 
to correct for the fact that we are including in the analysis 
AGN with logLx[erg s _1 ]<42.2 and wq to correct for this 
effect plus the redshift evolution of the AGN density: 



w(z) 



Wq(z) 



r- 



, Lx)dLx 



SI 
L 



,(*)■ 



(z,L x )dL x 



42.4 



0,L x )dL 



x 



oo 

Uir 



(z,L x )dL x 



(3) 



(4) 



<j)(z, L x ) is the soft XLF proposed in lEbrero etaLl (|2009f ) 
and Lu m depends on the survey flux limit and is a func- 
tion of redshift. I n detail, following E quations 9,11 and 
12 and Table 2 in lEbrero etaLl (|2009l ), we described the 
shape of the present-day luminosity function with slopes 
7i = 0.72 ± 0.02 and 72 = 2.04 ± 0.04, logL [h 7 oerg 
s _1 ] = 43.65 ± 0.05 which is the value of the luminos- 
ity where the change of slope occurs and normalization 
A = 3.76 ± 0.38 x 10- 6 lrf Mpc- 3 . We estimated the 
evolution factor assuming pi = 3.38 ± 0.09, P2 = —1.5, 
z c = 1.42, \ogL a = 44.6h 70 erg s" 1 and a = 0.100±0.005. 

Based on Eq. [3] and |4j the comoving space density rij 
of AGN hosting groups corrected for the XLF is given 
by: 



W J (Z,L X ) 



V, 



(5) 



max, j 



while the rij corrected for both the XLF and the z evo- 
lution is estimated by using: 
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Fig. 3. — Left Panel: Observed occupation of galaxy groups by XMM-COSMOS AGN as a function of the halo mass (open grey triangles) 
and after correcting for the soft XLF (green open circles) and for both the XLF and the AGN redshift evolution. Right Panel: Mean AGN 
occupation as a function of redshift for AGN hosting groups with M200 = 13 — 14Mq (colors with same meaning as left panel.) 



Fig. [5] (Right Panel) shows the mass function when all 
these effects are corrected. 

The ratio between the mass function of X-ray groups 
hosting AGN within -R200 by that of all X-ray galaxy 
groups generates the mean AGN HOD, which describes 
the occupation of DM halos by AGN. Fig. \5\(Left Panel) 
shows the observed average number of XMM-COSMOS 
AGN in a halo of given mass as a function of M200 (grey 
triangles) and the average number after correcting for the 
XLF (green open circles) and for both the XLF and the 
z evolution (black filled circles). The redshift evolution 
of the AGN fraction is shown in Fig. [2] (Right Panel). 
The grey circles represent the mean AGN occupation as 
a function of redshift for halos with M200 = 13 — 14M 
while the green circles and the black squares show the 
AGN fraction when we correct for the XLF only and for 
both the XLF and the z evolution effects, respectively. 
When we correct for the XLF, the mean occupation in- 
creases with z since in the AGN luminosity range of our 
interest, the AGN density increases with z, while the red- 
shift correction removes this trend producing a constant 
AGN fraction. 

We fitted the total AGN HOD assuming the model 
with a rolling-off power-law: 



We obtained constraints in the (a, f a ) - space and we 
found as best fit parameters of the mean AGN HOD, 
a = 0.06(-0.22;+0.45) and /„ = 0.06(0.04; 0.08), where 
the 68% confidence interval for a combined two parame- 
ter fit (Ax 2 = 2.3) is given in the brackets. Fig. \$\(Left 
Panel) shows the mean occupation of XMM-COSMOS 
AGN with the best fit parameters (solid black line) , and 
lcr confidence interval (shaded grey region), compared 
to the mean AGN HOD including C-COSMOS AGN 
in the analysis. The solid red line corresponds to the 
best fit model for the mean occupation of XMM+C- 
COSMOS AGN, with a = 0.06(-0.22, 0.36) and /„ = 
0.05(0.04,0.06), while the shaded region is the la confi- 
dence interval (see Fig. 21 Right Panel). 

Moreover by dividing the mass function of satellite 
(central) AGN host groups by that of all X-ray galaxy 
groups, we provide the fraction of AGN among satellite 
(central) galaxies as a function of halo mass (see Fig. [5j. 
We model the mean AGN occupation function in halos 
by decomposing it into the central and satellite contri- 
bution (N AGN )(M h ) = (N cen )(M h ) + (N sat )(M h ): 



i at mm \ f t ( logMh-logMn 

(N cen )(M h ) = f a erf 

\ &logM 



;-V. U ;.x)(-U,J=./„(^ 



exp 



M h 



(7) 



where f a is the normalization, M\ is the halo mass at 
which the number of central AGN is equal to that of 
satellite AGN, M cut is a cut-off mass scale. With our 
data alone, we cannot make meaningful constraints on 
Mi, ther efore we fixed logMi = 13.8M Q following the 
results of Miyaji et alJ f|201 ID - We verified that the re- 
sult does not change for logMi = 13 — 14.2M©, due to 
the fact that the fraction of AGN among central and 
satellite galaxies are comparable in this mass range. For 
the modelling of the rolling-off power-law, we assumed 
\ogM cut = 13.4Mq which corresponds to the mass below 
which our data points decay exponentially. 



(8) 



(N sat )(M h ) = f' a \ exp(-M cut /M h ) (9) 

where the central AGN occupation follows a softened 
step function and the satellite occupation a rolling-off 
power law (e.g. Kravtsov et a. 2004, Zheng et al. 2005, 
Zehavi et al. 2005, Tinker et al. 2005, Conroy et al. 
2006, Chatterjee et al. 2011, Richardson et al. 2012). In 
this formalism there are four free parameters f' a , M m i n , 
&iogM and a s , where M m i„ is the minimum mass where 
the occupation of central AGN is zero. 

As shown in Fig. [SJ the HOD of central AGN is de- 
scribed by a softened step function with logM m i„[M ] = 
12.75(12.10, 12.95) and a logM = 1.46(0.4,4.0) where the 
errors are the la confidence intervals estimated by us- 
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Fig. 4.— Left Panel: Occupation of galaxy groups by 41 XMM-COSMOS AGN (black open circles) and 58 XMM + C-COSMOS AGN 
(red filled circles) as a function of the halo mass, when correcting for the XLF and for the redshift evolution of the AGN density. The 
fit assuming a rolling-off power-law dependence of the HOD is shown as the solid black lines (best fit) and shaded regions (1<t confidence 
interval, A% 2 = 2.3). Right Panel: the confidence contours of the power-law best-fit parameters a and f a , for XMM-COSMOS AGN 
(black) and for XMM + C-COSMOS AGN (red) in galaxy groups. The contours mark the 68.3% and 95.4% confidence levels (respectively 
corresponding to A\ 2 = 2.3 and 6.17). 



ing the x 2 minimization technique with 2 free parame- 
ters. On the other hand, the satellite AGN HOD sug- 
gests a picture in which the average number of satel- 
lite AGN increases with M h (a s = 0.22(-0.07, 0.63) and 
f' a = 0.034(0.022, 0.046)) slower than the satellite HOD 
o f samples of g a laxies (cc M^ s=1 ~ 1,2 ). 

iMivaii et all (|2011ft first used a sample of ROSAT- 
RASS AGN and SDSS galaxies at z ~ 0.3 to study the 
occupancy of X-ray AGN in DM halos. They investigate 
three models: the first (model a) assumes that AGN only 
reside in satellite galaxies while the others explore the ef- 
fects of centrals (model b and c). In particular, in model 
b, the HOD of central AGN is constant and the satellite 
HOD has a power- law form at halo masses above M m i n . 
Model c has the same form except that only less massive 
DM halos contain central AGN. Our results clearly show 
that only a model which includes AGN in both satellite 
and central galaxies at any halo mass M > M m i„ can 
reproduce the observed AGN HOD. 

In agreement with our results, they found that the up- 
per limit of the power-law index of the satellite HOD is 
below unity, with a s < 0.95. In detail, their model b 
with best fit parameters a s — 0.55, logM min [h _1 M0] 
= 12.5 and f a — 2.9 x 10 ~ 2 is in perfect agreement 
with our results, showing that the luminosity and red- 
shift evolution of the mean AGN HOD is not strong in 
the luminosity and redshift ranges of our interest. In 
fact, we used an AGN sample with ( Lx) = 10 42 3 erg 
s _1 , while the results in IMivaii et al.l ([201 If ) provide the 
mean AGN HOD at (z) ~ 0.3 for more luminous AGN 
with (£0.1-2.4) ~ 10 44 ' 2 erg s _1 without any correction 
for the z evolution. Then the two models suggest a sim- 
ilar positive a s range, but negative values of the slope 
are not rejected (A^ 2 < 2.3) in their model. Note that 
while in our case the AGN fraction is a free parameter, 
they constrained it by normalizing the AGN HOD to the 



observed AGN number density. 

An assumption in the corrections applied by using Eq. 
[3] and IH is that the shape of the XLF remains the same 
between AGN in groups and the general AGN population 
at logL x [erg s _1 ]< 42.4. We still do not have sufficient 
observational basis to estimate the effects of possible dif- 
ference in the XLF shapes between group environment 
and the field. At higher luminosities, Krumpe et al. 
(2010; 2012) found that AGN at log(L x )[erg s" 1 ] « 44.6 
have higher bias values than those at » 43.9, correspond- 
ing to the mean host halo mass of log^MhHh ^MfT)] = 
13.0 and 13.3 respectively (jMivaii et al.ll201lj) . Viewing 
this trend from the Mh dependence of the XLF shape, 
the XLF in high Mh environments have XLF more biased 
towards higher Lx- If the positive correlations between 
Lx and the host halo mass Mh extends to luminosities 
in the range 41.5 < \ogLx [erg s _1 ] < 42.4, i.e. the 
range of logLn m (z) (see Fig. 1 left), the correction fac- 
tors using the overall XLF in Eq. 3 and 4 are overes- 
timated at higher Lx and underestimated at lower Lx- 
A quantitative assessment of this effect needs an esti- 
mate of the bivariate X-ray luminosity-host DMH mass 
function, which is far from being available. We com- 
ment that, if the higher Lx AGN preferentially occupy 
log M h [M ]> 13 DMHs than the field, our estimate of 
a s should be corrected to a lower value. 

4. TWO-POINTS STATISTICS 
4.1. Method 

The two-point auto-correlation function (ACF) de- 
scribes the excess probability over random of finding 
a pair with an object in the volume dV\ and another 
in the volume dV-2, separated by a distance r so that 
dP = n 2 [l + £,{r)]dV\dV2 : where n is the mean space den- 
sity. A known effect when measuring pair separations is 
that the peculiar velocities combined with the Hubble 
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Fig. 5. — Occupation of galaxy groups by satellite AGN (black 
squares) and central AGN (red crosses) as a function of the halo 
mass when correcting for the redshift evolution of the AGN density 
and for the soft XLF. The fit assuming a rolling-off power-law 
dependence is shown as solid black line (best fit) and dashed grey 
region (lcr confidence interval) for AGN among satellite galaxies. 
The HOD for central AGN has been modelled with a softened step 
function (solid red line) where the shaded orange region marks the 
68.3% confidence level. 

flow may cause a biased estimate of the distance when 
using the spectroscopic redshift. To avo id this effect it is 
usual ly computed the projected ACF (jDavis & Peebles 
1983): 



w p (r p ) = 2 



£(r p ,ir)dTr 



(10) 



where r p is the distance component perpendicular to 
the line of sight and ir parallel to the line of sight 



the line or sight < 
(|Fisher et alj fl994) 



The ACF is estimate d by using the m i nimum variance 
estimator described by lLandy fc Szalavl (|1993t ) : 



DD - 2DR + RR 
RR 



(11) 



where DD, DR and RR are the normalized number 
of data-data, data-random, and random-random source 
pairs, respectively. Eq. QT] indicates that an accurate es- 
timate of the distribution function of the random samples 
is crucial in order to obtain a reliable estimate of £(r p , it). 
Several observational biases must be taken into account 
when generating a random sample of AGN in a X-ray 
flux limited survey. In particular, in order to reproduce 
the selection function of the survey, one has to carefully 
reproduce the space and flux distributions of AGN, since 
the sensitivity in X-ray surveys is not homogeneous on 
the detector and therefore on the sky. An other impor- 
tant choice for obtaining a reliable estimate of w p (r p ), is 
to set iTmax in the calculation of the integral above. One 
should avoid values of ir ma x too large since they would 
add noise to the estimate of w p (r p ). If, instead, 7r max is 
too small one could not recover all the signal. 

We created an AGN random sample where each sim- 
ulated source is placed at a random position in the sky, 
with flux randomly extracted from the catalog of real 
source fluxes. The simulated source is kept in the random 



sample if its flux is above the sensitivity map value a t 
that position (jMivaii et al.ll2007t ICappelluti et al.ll2009D . 
The corresponding redshift for a random object is as- 
signed based on the smoothed redshift distribution of 
the AGN sample. 

Similarly, an unclustercd catalog of galaxy groups 
mimicking the selection function of the survey must be 
employed to quantify the degree to which the groups 
preferentially locate themselves in one another's neigh- 
borhood. We generate a random group catalog by calcu- 
lating at each area of a given sensitivity, the probability 
of observing a group of a given mass and redshift. We 
then used Monte-Carlo simulation of the group positional 
assignment, finally producing a catalog of hundred thou- 
sand objects. Also the X-ray surface brightness sensitiv- 
ity map is non-uniform in depth and consequently the 
probability of detecting groups of a particular mass is 
variable with redshift; in particular the minimum mass 
below which a group will be detected is an increasing 
function of z. 

4.2. Galaxy Groups ACF 

We measured the projected ACF of galaxy groups in 
the range r p = 0.1 — 40 h _1 Mpc by using Eq. [TU1 with 
Kmax = 80 h _1 Mpc (see FigfB]). The errors have been 
estimated using bootstrap resampling of the data, which 
consists of computing the variance of w p (r p ) in N rea i 
bootstrap realizations of the sample. Each realization is 
obtained by randomly selecting a subset of groups from 
the data sample allowing for repetitions. 

In the halo model approach, the clustering signal can 
be modelled as the sum of two contributions of pairs from 
the same DM halo (1-halo term) and those from different 
DM halos (2- halo term). In Fourier space, the 2-halo 
term can be explicitly written as (Seljak 2000, Cooray & 
Sheth 2002): 

P 2 - h ^b 2 P lln {k,z) (12) 



where Pu n (k,z) is the linear power spectrum and b is 
the bias factor of the sample. Then the galaxy groups 
two-point correlation function at large scales is given by: 



w 



V.2 



-h(r P ) 



b 2 

group 



2tt 



P Un {k)J (kr p )dk (13) 



where Jo{x) is the zeroth-order Bessel function of the 
first kind. Following this model, the galaxy group bias 
defines the relation between the 2-halo term of DM and 
groups clustering signal: 



Up,2 



-h{r p ) 



w DAi( r P, z = °) 



(14) 



According to this equation, we estimated an average 
bias factor in the range r p = 1 — 40 h Mpc equal to 
bgroup.obs — 2.20 ± 0.12 where the error corresponds to 
A\ 2 = 1 using a \ 2 minimization technique with 1 free 
paramete r. Following the bias m ass relation b[Mh, z) de - 
scribed in Ivan den Boschl (|2002t ) and lSheth etaLI (|2001l) . 
this observed bias, with respect to the DM distribution 
at z = corresponds to a typical mass \ogM typ [h _1 M0] 

- 13 fi [ S +0 07 

— lO.DO_ og . 

On the other hand the average bias of galaxy groups 
can be estimated starting from the known masses A/200- 
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Fig. 6. — : Upper Panel: Projected auto-correlation function of 
galaxy groups with z < 1. Following the halo model approach, the 
clustering signal at large scale (r p > 1 — 2h~ Mpc) is produced by 
separate DM halos (so called 2-halo term) and it can be expressed 



Fig. 7. 



Upper Panel: Projected cross-correlation function 



hi 



groups ]J ]\/J 



Op 



0) (dotted line), where the bias factor is 



defined in Equation 1151 Lower Panel: Bias factor versus r p . The 
solid line shows the best-fitting constant value. The shaded regions 
indicate the values of the bias for which A% 2 = 1 (dashed line) and 
Ax 2 = 4 (dotted lines). 

Since usually the bias factor is related to the average 
DM halo mass expressed in units of h~ 1 M@, hereafter 
we refer to Mh to indicate the galaxy groups masses in 
these units. Since we know the mass estimates, we can 
predict the bias factor of this sample of galaxy groups. 
By accounting for the fact that the linear regime of the 
structure formation is verified only at large scales, we 
estimated the average bias of the sample, including only 
pairs which contribute to the clustering signal at r p = 
1-40 Mpc h~ x . As described in lAllevato et al.l (|2011l ). 
we measured the average bias factor of the sample as: 



,(M h 



N, 



(15) 



pairs 



where b gg ^bggj is the bias factor of the group pair i — j 
and Np a i rs is the total number of group pairs in the 
range r p = 1 — 40 Mpc h~ x . The D factor is defined by 
Di(z)/Di(z = 0), where D-\ (z) i s the growth function 
(see eq. (10) in lEisenstein fc Hul (|2001| ) and references 
therein) and takes into account that the amplitude of the 
DM 2-halo term decreases with increasing redshift. 



group 



= 2.21 



+0.13 
-0.14 



By using this approach we obtained b 
where the errors have been estimated assuming a 20% er- 
ror on the galaxy groups masses. This value is in perfect 
agreement with the bias obtained from the ACF, with 

bgroup /bgroup,obs — 1.00 i 0.05. 

4.3. X-ray AGN-Galaxy Groups Cross- Correlation 

Measurements of the cross-correlation function be- 
tween AGN and groups use a vers ion of the estimator 
proposed bv lLandv fc Szalavl (|1993| ): 

1 



R„R A 



(D g D A - D g R A - D A R g + R g R A ) (16) 



of 253 XMM-COSMOS AGN with z s 



< 1 and 189 



iaxy 



groups, excluding from the analysis AGN that are within in galaxy 
groups. The clustering signal at large scale r p > 1 — 2h~ Mpc 
is due to AGN residing in different DM halos and is described by 

^ccf w dai( t v> z ) (dotted line), where the bias factor is de- 
fined in Equation [17] On the contrary the 1-halo term is zero since 
we are removing from the analysis AGN with i?200 responsible for 
the signal of AGN within the same DM halos. Lower Panel: Bias 
factor versus r v . The solid line shows the best-fitting constant 
value. The shaded regions indicate the values of the bias for which 
Ax 2 = 1 (dashed line) and A% 2 = 4 (dotted lines). 



where each data sample, with pair counts Di has an as- 
sociated random catalog, with pair counts Ri normalized 
by its number density. 

We estimated the cross-correlation function of 189 
galaxy groups and a subset of 253 XMM-COSMOS AGN 
with known spectroscopic redshift < 1, obtained exclud- 
ing those within i?200 ( see Fig. [7|). Following Eq. [T21 
the linear bias factor of the projected CCF of AGN with 
galaxy groups can be approximated by using the 2-halo 
term: 



bcCF,obs( r p) 



b A GN (rp)b group (r p 



w 



CCF 
p,2-h 



(r P ) 



w D M(r p ,z = 0) 

where b A GN and b grouv are the bias factor of AGN and 
galaxy groups, respectively and wdm is the projected 
dark matter CF. Fig. [7] (L ower Panel) shows the linear 
bias b 2 GGF as a function of r p over the scales r p ~ 1 — 
40 h _1 Mpc. We fitted the data points with a constant 
by using the x 2 minimization technique and we found 
b CCF obs — 3.90 ±0.28. The shaded regions show the bias 
values for which A% 2 =1 and 4 (68% and 99% confidence 
levels for one parameter). 



5. THE BIAS FACTOR IN THE HOD MODEL 

Based on the halo model approach, the AGN bias fac- 
tor depends on the mean AGN HOD (N AGN )(M h ): 



b A GN(z) = 



f£ . b h {M h ,z){N AGN ){M h )n h [M h ,z)dM h 



Jm ■ (N AGN )(M h )n h {M h7 z)dM h 



(18) 
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where M m i n is the minimum mass below which the AGN 
HOD is zero, nh(Mh) and bh( Mh) are the halo mas s func- 
tion and the halo bias given bv lSheth et aH(l200"l . Note 
that the large-scale bias factor does not depend on the 
normalization \ogMi[Mo\ — 10 138 and on f a . Since 
we are estimating the CCF excluding from the XMM- 
COSMOS sample 41 AGN in galaxy groups, the cluster- 
ing signal at large scale is due to AGN that live in halos 
with A/200 < 10 13 M Q or with masses that we can not 
observe at a given redshift z. This implies that the AGN 
bias can be written as: 



Jm* (z) b h (M h , z)(N AGN )(M h )n h (M h , z)dM h 



bAGN(z) 



J^(N AGN ){M h )n h (M h ,z)dM h 



(19) 

where M x {z) is the minimum mass we can observe for a 
group at redshift z (at z ~ 1 we only detect luminous 
and then massive groups). Note that we are assuming 
the separability of the mass and redshift dependence of 
the (N agn) j i.e. (N AGN )(M h ,z) = (N AGN )(M h ,z = 
0) x (N AG n)(z). 

Following this method, the bias factor b CCF is defined 
by: 

Ysi,jb A GN,ibgroup,jDiDj 



h 2 

°CCF 



where the sum is over the pairs i,j contributing to the 
clustering signal at large scale and b A GN(zi) is the AGN 
bias (Eq. HI?)) assuming a rolling-off power-law HOD 



based on our results. b group j is the bias a ssociated to 
the ga laxy group mass and redshift following Shet h et all 
(2001) and N pairs is the total number of AGN-group 
pairs in the range r p = 1—40 h _1 Mpc. We found 

bccF — 3.97lo]o5> which is in perfect agreement with 
the observed bias factor defined in Eq. [T7J The errors 
are due to the la errors on the power-law index a, M m i n 
and on b group . 

Similarly we can define the average halo mass (M) cor- 
responding to the observed CCF signal, i.e.: 



_ bAGN^bgroupjDjDjMl' 3 
/I f j bAGN,ibgroup,jDiDj 



(21) 



where M % h (z,i) is the average AGN host halo mass at red- 
shift Zi\ 



J^ {z) {N AGN )(M h )n h (M h ,z t )M h dM h 



M l h (z t ) 



lM^( N AGN)n h (M h , Zi )dM h 



(22) 



while Ml(zj) is the galaxy group mass, both in units 
of h _1 M Q . Fixing the values of a and M m i n , we found 
the average mass of XMM-COSMOS AGN to be log(M) 
[h- 1 M Q ]^ 13.20(13.10;13.25), where the errors corre- 
spond to the 68% confidence region in the \ogM m i n — a 
space. 

6. DISCUSSION 

Using a sample of X-ray selected AGN and galaxy 
groups in the COSMOS field at z < 1, we performed 
the first direct measurement of the AGN HOD in the 
mass range logM 2 oo = 13 — 14.5M , based on the mass 



function of galaxy groups hosting AGN. In contrast to 
previous works using the clustering signal of the sample, 
we directly counted the number of AGN within galaxy 
groups and we found 58 AGN in groups, associated to 22 
central and 36 satellite galaxies. This allowed us to put 
constrains on both the mean occupation of AGN among 
satellite and central galaxies as function of the halo mass, 
which p rovides informati o n on t he AGN triggering mech- 
anism. iStarikova et al.l (|201lD studied the halo occu- 
pation properties of AGN detected by the Chandra X- 
ray Observatory in the Bootes field over a redshift in- 
terval from z=0.17 to z=3, showing that X-ray AGNs 
are predominantly located at the centers of DMHs with 
M h > 4.1 x 10 12 h _1 Af Q , with an upper limit of the 
satellite fraction of 0.1 (Ax 2 < 2.3). The central loca- 
tions of the quasar host galaxies are expected in major 
merger models because mergers of equally sized galaxies 
preferentially occur at the c enters of DMHs (Hopkins et 
al. 2008). On the contrary iPadmanabhan et al.l (|2009f) 
observed the presence of the one-halo term in the cross- 
correlation function of optically selected z < 0.6 quasars 
and Luminous Red Galaxies, and use this to conclude 
that a large fraction of the AGNs is hosted by satellite 
galaxies. 

Our results show that the average number of AGN in 
satellite galaxies in the halo mass range ZogM/jM©] = 
13 - 14.5, and AGN luminosity log{L x ) [erg s _1 ]= 42.3 
might be comparable or even larger that the average 
number of AGN in central galaxies, i.e. X-ray AGN do 
not avoid satellite galaxies. A high fraction of AGN in 
satellite galaxies is expected in a picture where other 
phenomena like secular proc esses, might becom e dom i- 
nant in the AGN acti v ation. iMilosavlievic et all (12 006) . 
IHopkins fc Hernguistl ()2006| ) and IHopkins fc Hcrnguist 
(2009) showed that low-luminosity AGN could be trig- 
gered in more common nonmerger events, like stochastic 
encounters of the black holes and molecular clouds, tidal 
disruption or disk instability. This leads to the expec- 
tation of a characteristic transition to merger-induced 
fueling around the traditional quasar- Seyfert luminosity 
division. 

Moreover we found the power-law slope, which defines 
the evolution of the mean satellite HOD with halo mass, 
to be a s ~ — 0.6, suggesting a picture in which the 
average number of satellite AGN pe r halo increas e s with 
the halo mass. On the other hand, iMivaii et al.l ()2011f ) 
obtained a s < 0.95, but negative values of the slope are 
not rejected (Ax 2 < 2.3). 

It is interesting to compare this result with HOD analy- 
ses of galaxies. Previous HOD analyses of galaxies found 
a s ~ 1 — 1.2 for a wide range of a bsolute magnitudes 
and redshifts at least up to z ~ 1.2 (|Zehavi et al.l 120051 : 
Zheng ct al. 2007; Zcha vi et al.ll2010|) , implying a simple 
proportionality between halo mass and satellite number, 
(N sat ) oc Mh- Our results suggest that the mean HOD 
of satellite AGN might increase slower (a s < 0.63) with 
the halo mass respect to the linear proportion (a s = 1) 
in the satellite galaxy HOD, i.e. the AGN is not only 
triggered by the halo mass. On the contrary, a decreas- 
ing AGN fraction with the halo mass is consistent with 
previous observations that the AGN fraction is smaller 
in clusters than in groups in the nearby universe. 

In order to fully understand the growth history of 
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SMBHs as well as the physical processes responsible for 
the AGN activity we need to explore the AGN HOD at 
different redshifts, luminosities, and AGN types. The lu- 
minosity distribution of AGN that reside in halos of a 
given mass provides a tool to examine the distribution 
of halo mass for a given luminosity and study luminos- 
ity dependent clustering. While this formalism has been 
widely used in modelling galaxy clustering, it is still not 
applicable to AGN. In fact, it is also important to have 
larger numbers of AGN in galaxy groups which will en- 
able stronger constraints on the shape of the satellite 
and central AGN HOD, hopefully as function of AGN 
luminosity and redshift. 

7. CONCLUSIONS 

We have performed the first direct measurement of 
the mean halo occupation distribution of X-ray AGN 
as function of halo mass, by directly counting the num- 
ber of AGN within X-ray galaxy groups with masses 
logM 20 o[M Q ] = 13-14.5, in the COSMOS field at z < 1. 
Our findings are summarized as follows. 

1. We identified 41 XMM-COSMOS AGN within 
galaxy groups, defined as AGN located within 3 
times the group line-of-sight velocity dispersion 
and within -R200 and 17 additional sources includ- 
ing in the analysis C-COSMOS only selected AGN. 



2. We measured the mean AGN occupancy of galaxy 
groups as function of halo mass in the range 
logAf2oo[Afo] = 13—14.5 and we modelled the data 
points with a rolling-off power-law with the best fit 
index a — 0.06(0.22, 0.36) and normalization pa- 
rameter f a = 0.05(0.04, 0.06). 

3. Using a galaxy membership catalog, we associ- 
ated 22/58 and 36/58 AGN to central and satel- 
lites galaxies, respectively. We constrained that 
the mean AGN occupation function among cen- 
tral galaxies is described by a softened step func- 
tion above logM mm = 12.75(12.10, 12.95) while the 
satellite AGN HOD increases with the halo mass 
(a s < 0.63) slower than the satellite HOD of sam- 
ple of galaxies (a s = 1 — 1.2). 

4. We presented an estimate of the projected ACF of 
galaxy groups over the range r p = 0.1-40 h _1 Mpc 
at (z) — 0.5. We verified that the bias factor 
and the corresponding typical halo mass estimated 
from the observed galaxy group ACF, are in per- 
fect agreement with the values b group and (Mh) ob- 
tained by using the galaxy group mass estimates. 
In particular we found b group = 2.2lj^Q 14 and log 



(M h ) [rWg] = 13.61 



+0.09 
-0.10' 
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AGN IN GALAXY GROUPS 



ID 


Ra 




Dec 


z 


itJ & iw 200 


ID 


AON Ra 


AON Dec 


ID 


z c 


LogL/x 


Flap- 


Groups 


deg 




deg 


Groups 




XMM 


deg 


deg 


Gnandra 


AGN 


[h 70 erg s \ 


d 


11 


150.1898 


1 


.65725 


0.22 


14.28 


-99 


150.19864 


1.671924 


12390 


0.228* 


41.27 





17 


149.96413 


1 


.68033 


0.372 


14.02 


30744 


149.96396 


1.6805983 


-99 


0.372 


43.28 


1 


19 


150.37282 


1 


.60944 


0.103 


12.98 


2021 


150.37256 


1.6093965 


1678 


0.104 


41.22 


1 


20 


150.32494 


1 


.60313 


0.22 


13.38 


2186 


150.33601 


1.6012201 


1671 


0.234* 


41.88 





35 


150.2066 


1 


.82327 


0.53 


13.70 


-99 


150.20702 


1.823398 


1292 


0.529 


41.73 


1 


39 


149.82381 


1 


.8252701 


0.531 


13.63 


5502 


149.81248 


1.8238187 


219 


0.529 


42.37 





52 


150.44704 


1 


.8828501 


0.671 


13.64 


30681 


150.44731 


1.8832886 


690 


0.670* 


42.75 


1 


69 


150.42012 


1 


.9708 


0.862 


13.49 


5519 


150.42473 


1.9693011 


1185 


0.854* 


42.77 





78 


150.27461 


1 


.98884 


0.838 


13.61 


-99 


150.27823 


1.990302 


685 


0.838 


42.35 


1 


78 


150.27461 


1 


.9888 


0.838 


13.61 


206 


150.27434 


1.9884306 


1100 


0.673 


42.57 





79 


150.44728 


2 


.05392 


0.323 


13.77 


-99 


150.44765 


2.053961 


2605 


0.323 


41.62 


1 


87 


150.51109 


2 


.0269899 


0.899 


13.60 


2387 


150.51036 


2.0293686 


496 


0.899 


43.15 


1 


93 


149.6692 


2 


.0740 


0.338 


13.59 


417 


149.6692 


2.073962 


168 


0.34 


42.76 
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2 


.1103699 


0.361 


13.43 


6 


150.17978 


2.1101542 


42 


0.360 


43.19 


1 


118 


149.63463 


2 


.1357 


0.962 


13.70 


-111 


149.63828 


2.1494889 


323 


0.952 


43.25 
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2 


.2085 


0.186 


13.22 


302 


150.0571 
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1297 


0.186 


41.25 





124 


150.05656 


2 
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0.186 


13.22 


-99 
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1221 
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2 


.1811 
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2855 
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1634 
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2 


.15873 


0.377 


13.32 
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2.1551671 
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0.376 


41.24 


1 


130 


150.02382 


2 


.2032 


0.942 


13.84 


-99 


150.02303 


2.206567 


2413 


0.94 


42.45 
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2 
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0.676 
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42.08 
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41.44 


1 
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2 


.2614 


0.943 


13.97 
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149.50914 


2.261276 


-99 


0.962 


43.19 


1 
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150.28798 


2 


.2769 


0.122 


13.09 


317 


150.25331 


2.2779887 


3718 


0.165 


41.08 
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150.21454 


2 


.2801 


0.880 


13.77 


116 


150.20604 


2.2857771 


22 


0.874 


43.33 





149 


150.41566 


2 


.4302 


0.124 


13.85 


45 


150.33598 


2.4335926 


658 


0.121 


41.31 





149 


150.41566 


2 


.4302001 


0.124 


13.85 


73 


150.416 


2.4299896 


634 


0.124 


42.00 


1 


161 


149.95262 


2 


.3418 


0.941 


13.74 


430 


149.95949 


2.3560882 
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0.889 


42.85 
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149.66328 


2 


.2677 


0.676 


13.56 


5411 


149.66243 


2.2693584 


1646 


0.676 


42.76 





173 


150.05804 


2 


.38045 


0.347 


13.48 
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898 
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42.14 


1 
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149.63988 


2 


.3491 


0.950 
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149.64459 


2.3587193 
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0.962* 


42.88 
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2 
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13.55 
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41.72 
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2 
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2 
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2 
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o 
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2 
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2.642817 
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1 
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2 
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1 
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2 
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2 
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149.91588 


2.52139 


1310 
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1 
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2 
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60270 
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262 
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2 
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ollz 
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42.73 


1 
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2 
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13.57 
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1 
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1 
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24 
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150.0423 


2 


.6949 


0.219 


13.31 


5075 


150.04155 


2.6945302 


623 


0.221 


41.46 


1 


a Mass defined respect to 


200 times the 


mean density, with h 


=0.72. 















c = -99 means NOT detected. 

c Photometric or spectroscopic rcdshift. 

0: AGN among satellite galaxies; 1: AGN among central galaxies. 

Photometric rcdshift from Salvato ct al. (2011) . 



HOD of X-ray AGN 



13 



Sanders, D., Salvato, M., Aussel, H., et al. 2007, ApJS, 172, 86 
Schawinski, K„ Treister, E.,Urry, C. M., et al. 2011, ApJ, 727, 31 
Scoville, N., Abraham, R. G., Aussel, H., et al. 2007, ApJS, 172, 
38 

Shankar, F., Weinberg D. H., et al. 2009, ApJ, 690, 20 
Shankar F., et al., 2010, ApJ, 718, 231 
Shankar F., 2010, IAUS, 267, 248 

Shen Y., Strauss, M. A., Ross, N. P., Hall, P. B., et al. 2009, 

ApJ697, 1656 
Shen Y., 2009, ApJ, 704, 89 
Shen, Y., et al. 2010, ApJ, 719, 1693 
Sheth R. K., Mo H. J., Tormen G. 2001, MNRAS, 323, 1 
Silverman, J. D.; Kovac, K., Knobel, C, ApJ, 695, 171 
Silverman, J. D. et al. 2011, ApJ, 743, 2 
Starikova, S. et al., 2011, ApJ, 741, 15 

Tal, T., van Dokkum P. G., Nelan, J., et al. 2009, ApJ, 138, 1417 



Taniguchi, Y., Scoville, N. Z., Murayama, T., et al. 2007, ApJS, 
172, 9 

Tinker, J. L., Weinberg, D. H., Zheng, Z., Zehavi, I. 2005, ApJ, 
631, 41 

Ueda, Y., Akiyama, M., Ohta, K., & Miyaji, T., 2003, ApJ, 598, 
886 

Yang, Y., Mushotzky, R. F., Barger, A. J., & Cowie, L. L. 2006, 

ApJ, 645, 68 
van den Bosch, F. C. 2002, MNRAS, 331, 98 
Zamojski, M. A., Schiminovich, D., Rich, R. M., et al. 2007. 

ApJS, 172, 468 
Zehavi, I., et al. 2005, ApJ, 630, 1 
Zehavi, I., et al. 2010, arXiv: 1005.2413 
Zheng, Z., et al. 2005, ApJ, 633, 791 
Zheng, Z., Coil, A. L., & Zehavi, I. 2007, ApJ, 667, 760 
Zheng, Z.. Zehavi, I., Eisenstein, D. J., Weinberg, D. H., & Jing, 

Y. P. 2009, ApJ, 707, 554 



1 — 'Best Fit 



"i 1 1 r 



i 1 1 r 



"i r 



s o - 



Lowest a extreme 
Highest a extreme 



A 
•z 
u 
< 

V 
O 



-1 - 



-2 - 



-3 



11 




r 



J I U L 



J U L 



J I I L 



J L 



12 13 

log M h [h-i M sun ] 



14 



i 1 1 r 



c 

3 



3f2.5 



£ 12 

s 



4l.5 



log <M> = 12.5 



i 1 1 r 



X 



log <M> = 12.8 



J 1 1 L 



i 1 1 r 



log <M> = 13 



J I I L 



"i r 



-1 



-0.5 








1±J L 



a 



